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ABSTRACT 

We investigate a small-scale (=s 1.5 Mm along the slit), supersonic downflow of about 90 kms -1 in the transition region above the 
light-bridged sunspot umbra in AR 11836. The observations were obtained with the Interface Region Spectrograph (IRIS) on 2013 
September 2, from 16:40 to 17:59 UT. The downflow shows up as red-shifted “satellite” lines of the Siiv and Oiv transition region 
lines and is remarkably steady over the observing period of nearly 80 min. The downflow is not visible in the chromospheric lines, 
which only show an intensity enhancement at the location of the downflow. The density inferred from the line ratio of the red- 
shifted satellites of the Oiv lines (N t = lO la6±0 ~ 5 cm -3 ) is only a factor 2 smaller than the one inferred from the main components 
(N t = l0 10 - 95±0 - 20 cm -3 ). Consequently, this implies a substantial mass flux (« 5 x 10 -7 g cm -2 s -1 ), which would evacuate the overlying 
corona on time scales of the order of 10 s. We interpret these findings as evidence of a stationary termination shock of a supersonic 
siphon flow in a cool loop rooted in the central umbra of the spot. 
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1. Introduction 


Sunspots have been an area of intense research ever since Hale’s 
Hale| s ( 1908] ) discovery of strong magnetic fields in these struc¬ 
tures (see reviews by, e.g., |Borrero & Ichimoto|[2011| [Rem- 


pel & Schlichenmaier 201 1[ >. Here we investigate a particular 
sunspot phenomenon: a small-scale (~ 1.5 Mm along the slit), 
supersonic, relatively stable downflow of about 90 km s ' in the 
transition region above a sunspot umbra, which was first noticed 
by Tian et al. ( 2014a) > in the same data set. 

A similar phenomenon has first been seen in data obtained 
during the 21 July 1975 flight of the NRL High Resolution 
Telescope and Spectrograph (HRTS) by |Dere| ( 1982[ > and Nico¬ 
las et al. (1982i, and later in SOHO/SUMER observations by 
Brynildsen et al. ( 2001[ > and Brynildsen et al. ( |20Q4[ >. Since 
HRTS provided only a snapshot, no inferences could be made 
about the temporal variability or stability of this feature. Brynild¬ 
sen et al.j( |2001[ ), on the other hand, did have temporal informa¬ 
tion in the SUMER sit-and-stare time series they studied. In the 
20-min time series they presented they found line profiles that 
were well represented by two Gaussian line components in the 
N v 1242 A and O v 629 A lines, with the main component show¬ 
ing prominent 3-min oscillations and the red-shifted component 
showing nearly the same high velocity (w 90 kms -1 in one data 
set) without oscillations in both emission lines. They introduced 
the term “dual flow”. They did not discuss the origin of these 
dual flows in detail. Brynildsen et al. ( |2004 1 extended their ear¬ 
lier work and found “dual flows” (i.e. line profiles in which two 
distinct velocities could be observed within the same resolution 
element) in 5 out of the 12 sunspots they studied. 

There have been numerous reports of supersonic downflows 
above sunspots following the initial HRTS publications by Dere 
( ]1982| and |Nicolas et al.| ( | 1982] >. Most of them are based on 
observations with HRTS I and II and HRTS on Spacelab 2 


(e.g. 

Brekke et al. 1990 1991a b, 

Kjeldseth-Moe et al. 1988 

1993 

and often showed multiple flow components within the 


1" x 1" resolution element. It is not clear, however, whether 


all these flows are akin to the "dual flows" of |Brynildsen 
|et al. (2001). Gurman| ( |1993| l on the other hand measured only 
weak, mostly subsonic downflows when analysing observations 
of eight sunspots with UVSP on the Solar Maximum Mission 
(SMM). He noted an interesting aspect in which the sunspots 
of his work differed from the ones observed by HRTS: the latter 
ones all had obvious light-bridges, while only one of the eight he 
observed with UVSP did. |Kjeldseth-Moe et al.| ( [T988[ ) did also 
speculate that the flows, which they found to have a typical size 
of 5", may be associated with structures such as light bridges. 

In addition to the two Brynildsen et al. papers cited above, 
there have been several other studies of transition region down¬ 
flows above sunspots based on data from SUMER and CDS on 
SOHO (e.g. |Brynildsen et al.||1998[ |Brosius & Landi||2005[ l. 
These downflows, however, we believe are different from the 
phenomenon we are describing in this paper. They were not 
as strong, mostly subsonic or in any case below 35 kms ', and 
more extended, often covering a good fraction of the penumbra. 


Supersonic downflows in the transition region above a 
sunspot have also been recently detected by |Kleint et al.| (|2014|> 
with the Interface Region Spectrograph (IRIS; De Pontieu 


|et al.||20f4) ). These authors observed bursts of extremely broad 
line profiles in Mgn h and k, Cn 1336A, Siiv 1394A and 
1403 A, suggesting supersonic downflows of up to 200kms 1 
and weaker upflows that are associated with small-scale bright- 
enings in sunspot umbrae and penumbrae. The broad line pro¬ 
files do not have distinct velocity components as e.g. the multiple 


flows observed by 

Dere |1982|, Nicolas et al. (1982 ), 

Kjeldseth- 

Moe et al. (1993), 

Brynildsen et al. (2001), and Bryni 

dsen et al. 

(2004) and they appear rather intermittently in short 

lursts last- 
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ing about 20 s. These characteristics are very different from the 
rather steady, “mono-chromatic” downflow discussed in this pa¬ 
per, which reveals itself as Gaussian “satellite” lines in the spec¬ 
tra of Si iv and Oiv, well separated from the main components. 
Interestingly, such a well-separated component is also visible in 
the T-t-plots displayed in Figure 3 of Kleint et al. (2014 1 at 
about 100 km s 1 red-shift, but not discussed in the paper. 

In Section 2 we will describe the observations. The analysis 
and results are presented in Section 3 and discussed in Section 
4, before concluding in Section 5. 


2. Observations 


For this study we used the same data set as was used by Tian et al. 
( 2014a] > for their study of the propagation behaviour of shock 
waves in a sunspot. The data were obtained with the IRIS in sit- 
and-stare mode on 2013 September 2 from 16:39:35 to 17:58:45 
UT. The target was the sunspot in AR 11836 near disk center at 
x=106", y=58"(coordinates of the center of the slit at the be¬ 
ginning of the observations). Slit length and width were 68"and 
07166, respectively. The pixel size along the slit is 07166. The 
data were acquired by binning the FUV data in the wavelength 
direction by a factor of 2, and thus have a dispersion of ~ 25 mA 
per pixel (about 5 km s _1 at 1400 A) in both the NUV and FUV 
channels. The cadence of the spectral data was 3 s, with expo¬ 
sure times of 2 s. The observations comprise nine spectral win¬ 
dows: Cn 1336 A, Fexii 1349 A, Cli 1352 A, Oi 1356 A, Siiv 
1394 A, Si iv 1403 A, NUV at 2786 A and 2831 A, and Mgn h 
and k 2796 A. The spectral window containing the Si iv 1403 A 
line also includes the Oiv 1400A, 1401 A, and 1405 A lines, 
the last one blended with a S iv line. Parallel to the spectral data 
slit-jaw images were taken in three filters (2796 A, 1400 A and 
1330 A) at a 12 s cadence. The pixel size of the slit-jaw images 
is 07166x07166. We used level-2 data from the IRIS archive at 
Lockheed-Martin Solar and Astrophysics Laboratory (LMSAL), 
i.e. data which have already been corrected for dark current, flat 
field and geometrical distortions. Below we also discuss the line 
ratios of the Si iv doublet and the O iv density sensitive lines. For 
that purpose we applied the current estimate of the wavelength- 
dependent effective area as provided by the Solarsoft routine 
iris_get_response. pro, using version 3 of the calibration, 
released on 2015 April 1. 

For context and to gain some information about the higher 
temperature plasma and the magnetic field structure we comple¬ 
mented the IRIS observations with data from the Atmospheric 
Imaging Assembly (AIA; Lemen et al .|2012 1 and the Helioseis- 
mic and Magnetic Imager (HMI; |Scherrer et al.||2012|l on the 
Solar Dynamics Observatory (SDO; Pesnell et al.|2012 l. To al¬ 
low a comparison with the IRIS observations, the AIA and HMI 
data were remapped to the IRIS data, both spatially and tempo¬ 
rally, using metadata in the FITS headers and the map routines 
of SolarSoft. No additional roll angle correction has been ap¬ 
plied. A secondary step in SDO-IRIS cross-alignment has been 
achieved by comparing the HMI white-light images and the near 
continuum data from the 2831 A spectral window of IRIS (wave¬ 
length index 15). The different IRIS spectral windows have been 
cross-aligned using the fiducial lines. No further inter-instrument 
alignment between HMI and the different AIA channels has been 
applied. In a final step, temporal jitter and drift of IRIS along the 
slit has been removed from the spectral data by cross-correlation 
techniques applied to the near continuum intensity in the 2831 A 
spectral window. The totally accrued drift during the observa¬ 
tions was less than 3.5 IRIS pixels or, equivalently, one HMI 
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pixel (075). Pointing excursions of IRIS perpendicular to the slit 
cannot be corrected for. Altogether, we believe that the cross¬ 
alignment between the IRIS and SDO data is within one HMI 
pixel. 


3. Observational results 

We start with a summary of the characteristics of the observed 
supersonic downflow. The various findings are discussed in more 
detail below. Note that throughout this paper we adhere to 
the convention that positive velocities represent blueshifts. The 
downflow 

- shows up as red-shifted “satellite” lines of the Siiv and Oiv 
transition region lines, completely separated from the main 
components by about 90km s -1 (Fig. [T|—14]), 

- is very localised and small-scale (~ 1.5 Mm along the slit; 

Fig# 

- occurs in the center of a sunspot umbra, which has one 
prominent and two smaller light-bridges (Fig. [Tji, 

- is a very steady-state flow both in line intensity and veloc¬ 
ity (Fig. [4 1 , lasting at least as long as the observations, with 
surprising y little variation for the transition region, which 
usually is highly dynamic (cf. Kleint et al., 2014), 

- does not participate in the 3-min shock wave dynamics, 
which dominate the main components (Fig. [4]), and shows 
only very weak, anti-correlated intensity signatures (Fig. [6|, 

- is comparatively more pronounced (at about 60% ± 20% of 
the main component) in the O iv lines than in the Si iv doublet 
(~ 20% ± 10% of the main component), 

- has a density of about 4xl0 ll) cm~ 3 , with little temporal and 

spatial variations, and accordingly a mass flux of ~ 5 x 10 7 g 

-2 -1 

cm s , 

- has a density about a factor ~ 2 lower than the one inferred 
from the main component (Fig#, 

- is visible only in the Siiv and Oiv transition region lines, 
with no velocity signature in the chromospheric lines Mg n h 
and k or Cii (Fig#, 

- is co-spatial with an increased intensity in those same 
chromospheric lines, suggesting related enhanced heating 
(Fig#, 

- is not characterised by distinctly peculiar line profiles; in par¬ 
ticular the profiles of the Siiv and Oiv satellite lines are 
nearly Gaussian with widths which are about the same as 
that of the main components, although the width of the Si iv 
satellite lines tends to be larger than the main component, by 
about 40% with large uncertainties (Fig#. 

3.1. Context of observation 

Figure |T] provides an overview of the region studied in this pa¬ 
per. The upper left shows a HMI full disk magnetogram, and 
the upper right a cut-out HMI magnetogram and AIA 171 A in¬ 
tensity image of AR 11836. The lower left shows an even more 
zoomed-in HMI intensity image of the central sunspot of AR 
11836. The position of the IRIS slit is indicated (x = 0). The 
lower right shows a 1-min average of the spectral window con¬ 
taining the Si iv 1403 A line. The red marks indicate the region 
(about 2.8 Mm wide) in the umbra where the Si iv and O iv lines 
show satellites. In the northern part the satellites are strongest; 
the size of this northern part is about 1.5 Mm. The center of this 
region is labelled “C” and used as the reference point of the satel¬ 
lite feature in the following. 




















Straus et al.: A Steady-state Supersonic Downflow in the TR above a Sunspot Umbra 




-20 -10 0 10 20 
rel. x [Mm] 


1400 1401 1402 1403 1404 1405 

Wavelength [A] 


Fig. 1 . Upper panels (context images of AR 11836): full disk (left) and cut-out (middle) HMI magnetogram and corresponding AIA 171 A image 
(right). Lower panels: zoom of HMI continuum image on the area marked by the red box in the upper panel, with the IRIS slit marked at x=0 
(left), and the corresponding spectra in the Siiv 1403 A window averaged over 1 minute (right panel). All data were taken on 2013 Sep 2, 17:14 
UT. Regions where the slit crosses the penumbra (PU) and umbra (U) are marked. The central part of the umbra where the Si iv and O iv lines 
show satellites is marked in red. In the northern part the satellites are strongest, the center of this region is labelled “C” and used as the reference 
point of the satellite feature in the following. The prominent upflow labeled "T" is due to a short-lived transient. 


The same data set has been studied and described by |Tian| 
et al. (2014a|l and Yurchyshyn et al. (2015 1 , who investigated the 


strong 3-min umbral oscillations found in this spot. It is interest 
ing to note that point “C” seems to coincide both with the appar¬ 
ent source of the umbral running waves studied by Yurchyshyn 


et al. (2015) and with the location of the most intense magnetic 


held in the sunspot. 


As can be seen in the HMI images, the spot of AR 11836 is 
an isolated, small, round spot of leading polarity. There is no fol¬ 
lowing spot, only active region plage. The spot is surrounded by 
a large area with network field of the same polarity. Opposite po¬ 
larity can be found in onl y a few locations at the south and north¬ 
east of the spot. The loops visible in AIA 171 A are inclined and 
connect the spot to plage of opposite polarity following the spot. 
The loops seem to be anchored in the dark north-west quadrant, 
whereas the rest of the sunspot displays a bright plume in the 
AIA 171 A images. 


The sunspot shows an irregular umbra, split up in two pieces 
by a prominent light-bridge in the western part. The bigger, east¬ 
ern part is further split into three pieces by two smaller, less 
prominent light-bridges. According to Gurman ( 1993) > it might 
be an important aspect for the presence of a supersonic down¬ 


flow that this sunspot had a light-bridge. The slit of the IRIS 
spectrograph crossed the eastern part centrally in the north-south 
direction. 


3.2. Average properties 

The Si iv and Oiv lines form at approximately 70 000 K and 
150 000 K, respectively. In both lines, the supersonic down¬ 
flow is stable enough to be prominent in the spectra averaged 
over the whole 80 minutes of the observations (see Fig.[2|. The 
strong northern satellite region is the darkest part in the “pseudo¬ 
continuum” between the Mgn h and k lines and shows at the 
same time the brightest intensity in all emission lines. There is no 
hint of any red-shifted components in the chromospheric spec¬ 
tral bands, in particular also not in the strong Mg n h and k and 
C ii 1336 A lines. 

In order to proceed with a quantitative analysis of the line 
profiles, we determined the main line parameters, namely the 
peak value, 7 max , its position, I rnax , and the full width at half max¬ 
imum, FWHM, by means of a parabolic fit to the central three 
pixels of each profile. Since the core of the Siiv and Oiv lines 
is only 3 or 4 pixels wide in this data set, a fit to a larger num- 
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Fig. 2. Spectra averaged over the whole time series as function of wavelength and position along the slit (lower panels, logarithmic color scale) 
and further averaged over the different regions (upper panels). The marks of different regions in the lower row correspond to those given in Fig. [T] 


ber of pixels is normally not required. At the extrema of shock 
maneuvres, double-peaked profiles may appear: the fit procedure 
then identifies the strongest peak. If needed, the total line inten¬ 
sity, / to t, is computed like in the case of Gaussian profiles as the 
peak value times the line width: / t01 = 0.5 VO 71 '/ log 2) x 7 max x 
FWHM = 1.064467 x 7 max x FWHM. We have compared this 
procedure with a more standard Gaussian fit, verifying that the 
agreement is excellent, especially for the more prominent lines. 
The Gaussian fit procedure, however, tends to fail more often in 
the case of weak profiles. The double peaks at the extrema of 
shock maneuvres are also less clearly identified. Finally, Gaus¬ 
sian fits assume symmetric profiles, while parabolic fits do not 
make this assumption. 


The average properties of the line parameters thus estimated 
were then determined from the histograms computed over 13 
pixels centered on point “C” along the slit, i.e. the region within 
the two marks nearest to point “C” in Fig. |T](the northern satel¬ 
lite region). A Gaussian fit to the peak of the histogram yields 
both the average value and the standard deviation. We estimated 
the averages both over the full time span and for the first 30 min¬ 
utes of the observing run only. We computed the average over 
the first 30 minutes because after that time the red-shifted com¬ 
ponents become fainter. Moreover, after about 30 minutes, the 
mean position along the slit of the downdraft region starts drift¬ 
ing south with respect to point “C” by several pixels. The average 
values do not change much in either case: for the Si iv line the 
mean shift between the two components is 91 + 8 km s' 1 , while 
the average value for the first 30 minutes is 94 + 6 km s 1 ; the 
corresponding values for the Oiv lines are 82 + 8 km s' 1 , while 
the average for the first 30 minutes is 84 + 6 km s 1 . In the fol¬ 
lowing we will always refer to averages computed for the first 
30 minutes, unless explicitly stated otherwise. 
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The satellite lines are most pronounced in O iv, with the total 
intensity being on average approximately 60+15% of the main 
component. The average Siiv satellite line is about 20+10% of 
the total intensity of the main component. Variations in the rela¬ 
tive intensity of the satellite to the main component in different 
O iv lines can be interpreted as due to density variations (see dis¬ 
cussion of Sec. 3.3.2) . Since we do not have estimates of plasma 
density in the region where the Si iv lines form, we are unable to 
offer a simple explanation of the different behaviour of this ra¬ 
tio between satellite and main component to the Oiv lines. This 
piece of information may be useful as constraints to more de¬ 
tailed models of this phenomenon. However, we note that al¬ 
ready Nicolas et al. (1982), using a larger set of lines, estimated 
that the shape of the differential emission measure of the super¬ 
sonic downflows is drastically different from the curve corre¬ 
sponding to the surrounding umbra. In particular, the differential 
emission measure they found in the downflows (their Fig. 8) is 
similar to the one they find in the sunspot umbra for log T > 5, 
while it continues to drop towards lower temperatures, as op¬ 
posed to the increase after that turn-around point seen in the 
umbra. This implies that the intensity of the O iv in downflows 
should be comparable to the one observed in the surrounding 
umbra, while the Si iv lines intensities should be strongly sup¬ 
pressed, as observed. 


We performed a similar analysis on the southern part of the 
satellite region. There the satellite lines are weaker compared to 
the main component: the Oiv satellite lines are on average only 
about 10% of the main component, while the Siiv satellites are 
5% or less. However, the velocity shifts remain as stable and 
comparable in magnitude to the values measured around point 
“C”, although slightly smaller on average: ~ 70 km s for the 
O iv lines, and * 80 km s 1 for the Si iv lines. 
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OIV OIV Si IV OIV/SIV 



Fig. 3. Symbols: 1-min average line profiles 
at point “C” with error bars (standard deviation 
of the twenty individual measurements), taken 
5 minutes before the reference spectrum in Fig. 
[T]( 17:09 UT). At this time, the satellite lines in 
Oiv are particularly strong. Solid line: average 
spectrum outside the spot. The main and satel¬ 
lite lines have all similar line widths (FWHM=s 
30kms _1 ), with the exception of Sirv 1403 A 
which has a significantly narrower main com¬ 
ponent. The satellites are less pronounced in 
Si iv than in O iv. 


3.3. Plasma diagnostics 

For the remainder of this section, except when noted otherwise, 
we analysed line profiles calibrated according to the procedure 
described in IRIS Technical Notes #1 and #244 and adopted the 
most recent effective areas for the IRIS spectrograph as provided 
by the SolarSoft routine iris_get_response .pro (version 3, 
2015 April 1). 


3.3.1. Line optical thickness 


In the plage regions outside the sunspot, the profiles of the C ii 
lines at 1334.5 A and 1335.7 A display an evident self-reversal, 
a clear signature of a significant optical thickness. Even in the 
penumbra the ratio of the total intensities of the two components, 
in contrast with the optical thin limit of 1:2, is around 1:1.4, a 
value close to the observed value in quiescent regions ( (Judge| 
et al. 2003) and to the computed values in the reference quiescent 


chromosphere by Avrett & Loeser (2008!. 

On the other hand, the profiles in the region around point 
“C” are nearly Gaussian. Nevertheless, the line widths of the C ii 
1334.5 A are on average larger than those of the Cn 1335.7 A 
line in the northern satellite region (1.3 ± 0.2) while at the same 
time the ratio of the total line intensities is consistent with the 
optical thin value (2.0 ±0.3^] These two measurements together 
suggest that the C n doublet in that region is affected by a modest 
optical thickness which is possibly smaller than in the quiet Sun, 
and that it is most likely being formed in the effectively optically 
thin regime: the line profiles become broader before the total line 
intensity is reduced by the destruction of photons due to colli¬ 
sions in the multiple scattering process.These measurements are 
consistent with earlier estimates (e.g. Jordan et al.|1978 Bartoe 


et al. || 1979) Tian et al. 2009) suggesting a lower opacity of the 
chromosphere above sunspot umbrae. 


1 http://iris.lmsal.com/documents.html 

2 The 1335.7 A line is in reality composed of two lines, at 1335.71 A 
and 1335.66 A; this latter line, however, is much weaker and does not 
therefore affect significantly the width of the observed profile. 


The average line profiles of the Siiv 1393.8 A and 1402.8 A 
lines do not show any sign of self-reversal or of significant de¬ 
partures from Gaussian in the line core in any region crossed by 
the slit (Fig. [2]»; we could not find any clearly non-Gaussian pro¬ 
files that could be attributed to optical thickness effects in the 
data we browsed, although we could not check all the individual 
profiles in the data set. 


A first analysis of the intensity ratio in the northern satellite 
region around point “C” of the main components of the two Si iv 
lines showed nevertheless the mean ratio of the 1394 A to the 
1403 A to be around 1.4, a value which is significantly lower 
than the optical thin line ratio of 2. During the course of the revi¬ 
sion of this paper, however, an updated, post-launch radiometric 
calibration of IRIS was released. With this latest calibration, the 
ratio of the main components of the Siiv 1393.8 and 1402.8 A 
lines became 1.9+0.1, much closer to the optically thin ratio. The 
ratio of the widths of the main component of the Si iv 1402.8 A 
line to the 1393.8 A line is on average 1.04 ± 0.08 
region, a value close to the value of 1.07 given by 
|Feldman| ( |2004| l for the quiet Sun. 


in the same 
Doschek & 


If we regard the results of the analysis of the C n lines above 
as relevant also for the main Si iv lines, i.e. that the optical thick¬ 
ness in point “C” is even smaller than in the average quiet Sun, 
we can also cite previous studies carried out in quiescent re¬ 
gions. In particular, the analysis of Roussel-Dupre et al. (1979) 
on the limb brightening of the line intensity of Si iv 1393.8 A in 
OSO-8 data indicates an optical depth of 0.22 ±0.16 in the quiet 
Sun, in contrast with the much larger optical thickness of the C n 
1334.5 A and 1335.7 A doublet estimated by Avrett & Foeser 


( 2008) 1 . Moreover, the Siiv lines 1393.8A and 1402.8A lines 
are not listed in Table 3 of Doschek & Feldman (2004 ) among 
the lines with measurable optical thickness in the FUV spectral 
range observed by SOHO/SUMER. 


On the other hand, the ratio of intensities of the satellite lines 
with the same post-launch calibration has the lower value of 
1.6 ± 0.3. The ratio of widths is 1.1 ± 0.3. Moreover, the Siiv 
satellite lines are clearly broader than their main components: 
the ratio of the FWHM of the satellite lines to the corresponding 
main component is 1.35 + 0.4 (see also Fig. |3j. This is in contrast 
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Fig. 4. Upper left panel: temporal average of AIA 304 A intensity, with HMI continuum intensity contours from Fig.[I]overlayed. The green solid 
and dashed lines indicate a path of inflows and the area of integration used for the time-distance plot in the right panel. Lower left panels: temporal 
average images of HMI continuum intensity and three AIA channels. The red marks indicate the satellite region and point “C” as in Fig. [T] Central 
panels: temporal evolution of the Oiv and Siiv lines averaged over 13 pixels in the northern satellite region centered on point “C”. The region 
of the satellite in Siiv, marked by the bar on top, has been contrast-enhanced to better show the weak connection to the main component’s shock 
maneuvres. Right panel: relative AIA 304 A intensity fluctuations integrated across the path marked by the green lines in the context image in the 
upper left. The colortable maps variations from -30% to +30%. Coronal rain is visible with plane-of-sky velocities of the order of the downflow 
velocity of the satellite (90 kms -1 , marked by the white line). However, these inflows do not seem to extend to the position of the satellite feature 
in the umbra, but appear to end at the light-bridge as marked by the vertical, blue lines in the time-distance plot and the corresponding marks on 
the path in the HMI continuum context image. 


with the observation that the widths of the O iv satellite lines are 
approximately equal to those of their main components (ratio: 
0.95 + 0.15). Since the peak intensities of the Si iv and Oiv satel¬ 
lite lines are comparable (see again Fig. [3], it is unlikely that a 
bias due to a possible instrumental background or solar contin¬ 
uum could be invoked to explain these measurements. 


It would then be tempting to interpret the larger width of the 
Si iv satellite lines as an optically thickness effect. Following the 
approach of |Doschek & Feldman| ( [2()04j i, this increased width 
would imply an optical thickness of the order of unity or larger 
(see for instance their Fig. 3). An optical thickness of the order 
of unity in the Siiv 1394 A line, with the densitie s of th e order 
of 10 1 * °— 10 11 cm~ 3 as inferred in the following Sec. 3.3.2 implies 
a geometrical thickness of the emitting plasma of at least a few 
hundreds of kilometers. Whether optical thickness effects alone 
can explain the broader Si iv satellite lines (as well as their ratio) 
will need to be verified by more detailed models. 


3.3.2. Electron density 


The intercombination multiplet of Oiv lines at 1397.20A, 
1399.77 A, 1401.16 A, 1404.81 A, and 1407.39 A provides a 


well k nown set of density-sensitive pairs (e.g.:|Dwivedi & Gupta 


1992); the lines at 1397.20 A and 1407.39 A are however outside 


the spectral window of this data set, while the line at 1404.81 A 
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is blended with a S iv line. We thus considered the ratio between 
the Oiv 1399.77 A and 1401.16 A lines as density diagnostics. 
We adopted the CHIANTI database of atomic data version 7.1.4 
( |Dere et ak 1997| Landi et al. 201 3[), assuming a Maxwellian 

distribution. It should be noted that |Dudrk et al. ( 201 4[ ) recently 
studied the effects of non-Maxwellian distributions of electron 
energies on both the Siiv and Oiv multiplets. We feel however 
that it would be unlikely that the stationary flows we are examin¬ 
ing could sustain substantial departures from equilibrium energy 
distributions. 

The time variation of densities inferred from this ratio is 
shown in Fig. [5] alongside the variation of the line intensities 
of both the main and the satellite components. To improve the 
signal-to-noise ratio, the line intensities were smoothed using a 
bi-Gaussian kernel with FWHM (full width at half maximum) of 

3 pixels along the slit and of 5 spectra in the time domain, corre¬ 
sponding to 0.35 Mm and 15 s, respectively. Note that smoothing 
data does not resolve the issue of the presence of possible errors 
in determining the background continuum which is to be sub¬ 
tracted from the line intensities before computing the line ratio. 
Since the Oiv 1400A line intensity is on average about 1/4 of 
the O iv 1401 A line in the region showing the supersonic down¬ 
flow, a systematic error in the background correction will affect 
the former line more than the latter. The result would be a sys¬ 
tematic bias in the determination of the density. For instance, 
a residual background intensity in the O iv lines would explain 
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Fig. 5. Variation of the intensity of the main component of the Oiv 1401 A line in the region that showed a supersonic downflow (first panel from 
the top), together with density inferred from the ratio of the Oiv 1400 A and 1401 A lines (second panel). The variation of the intensity of the 
satellite line at 1401A (third panel) and the density from the ratio of the satellite lines at 1400 A and 1401 A are shown in the third and fourth 
panel, respectively. The horizontal lines mark the boundaries of the regions given in Fig.[T] the dashed line marks point “C”. 


why the line ratio and therefore the density appears to increase 
in those regions where the lines are weaker (e.g. at the edge of 
the region where the satellite lines are visible). 

We examined the statistical distribution of densities derived 
from the ratio of the main components of the two O iv lines for 
the full spatial and temporal range over which the downflow is 
observable in this data set, i.e. for the full 80 minutes of the ob¬ 
serving run and within 2.5 Mm from point “C” (Fig. |5J. We 
obtained a mean value of log N e = 10.95 ± 0.20 (1 <x uncer¬ 
tainty, N e in cm -3 ), while the densities from the satellite lines 
are \ogN e = 10.6 ± 0.25. 

3.4. Signatures of shock dynamics 

The temporal fluctuations of the Siiv and Oiv lines show their 
main components with prominent shock wave passing maneu- 
vres and rather stationary supersonic downflow satellite lines. 
The satellite lines are not participating in the shock wave ma- 
neuvres, but in both lines the satellite intensity is slightly 
linked to the Doppler shift of the main component, showing 
a slightly brighter satellite feature when the main component 
passes through zero velocity in the steep rise of the shock front 
from maximum red-shift to blue-shift (Fig. [4] and [6}. This sug¬ 


gests that the plasma of the downflow is to a first approxima¬ 
tion dynamically decoupled from the plasma that forms the main 
components. One could speculate about a plume of transition 
region temperature plasma above the actual transition region. 
This would be a different scenario from the cartoon presented in 
Fig. 14 of Nicolas et al. ( |1982) >. However, we detect subtle sig¬ 
natures of a connection between the main and red-shifted com¬ 
ponent, as we discuss below. 

In Fig. [6| the main line parameters determined as described 
in Sec. |3.2] are shown after averaging over all 13 pixels in the 
northern part of the downflow region. The shock-like temporal 
fluctuations are striking in both the Doppler shift (black curve in 
upper right panel) and main peak intensity (red curve), especially 
in the middle of the displayed part of the observations. The green 
crosses mark the zero-crossing times of velocity in the shock 
front passages. At those times the main peak intensity is near its 
minimum. 

The symbols in the lower right panel of Fig. [6] show the in¬ 
tensity of the Si iv 1403 A red-shifted satellite as a function of 
time. Here we averaged the satellite peak intensities over five 
consecutive measurements in order to reduce the noise. One can 
see that the shock waves leave a weak imprint on the satellite 
intensity. However, there is an anti-correlation to the intensity of 
the main component. The satellite’s maximum intensity is found 
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Fig. 6. Right panels: temporal evolution of Si iv line parameters averaged over the region where the line shows a satellite (black: intensity; red: 
velocity). Dotted and dashed vertical lines mark the beginning and end of shock front passages, respectively. Positive velocities correspond to blue 
shifts. The times of zero velocity in the shock fronts are marked by green crosses. The main component’s peak intensity shows a minimum at 
these times. The lower right panel shows the temporal evolution of the satellite’s peak intensity. Left panels: temporal evolution of the average 
shock front of the 9 shocks marked by crosses in the upper right panel; time is given relative to the zero-crossing times of the main component’s 
Doppler-shift. The intensities of the Siiv satellite and main component are anti-correlated. 


at times of mimimum intensity in the main component (see black 
symbols in the left panels of Fig. [6] which show the fluctuations 
averaged over the 9 shock front events having their zero-velocity 
marked by crosses in the upper right panel, taking the time of 
vanishing velocity as temporal reference). The error bars in the 
left panels show the error of the mean value of the displayed 
signal after averaging over the 9 shock events. 

Interestingly, the satellite intensity (black symbols in the 
lower left panel of Fig. [6]» is not only anti-correlated to the in¬ 
tensity of the main component, but also shows a different shape. 
The shock-like shape is much less pronounced in the satellite’s 
intensity profile and reveals a slight preference of opposite sign 
having a slow rise before and a faster decrease after the peak in¬ 
tensity, which might hint at a downward shock. We will come 
back to this particular feature further down in the Discussion. 


4. Discussion 


The transition region is highly structured and dynamic (e.g. 
Mariska 1992). This fact has become even more obvious in re¬ 
cent high-resolution IRIS observations (e.g. Tian et al. [2014b 


Hansteen et aTj j2014| . It is therefore quite remarkable, that the 
downflow we observe in a sunspot umbra is rather stationary for 
nearly 1.5 hours. Has such a phenomenon been observed be¬ 
fore? We do believe that the one snapshot obtained during the 
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HRTS I flight by Dere ( 1982] ) and Nicolas et al. ( |1982 1 and the 
SUMER observations of a dual flows in O v and N v in sunspots 
by |Brynildsen et al.] ( |2001[ ) and |Brynildsen et ak] ( |2004| ) show the 
same phenomenon, albeit at a different (higher) temperature in 
the latter cases (O v and N v vs Si iv and O iv). 


What could be the origin of the feature we see? An obvi¬ 
ous candidate, of course, is coronal rain, a commonly observed 
phenomenon, in the EUV and even Her (e.g. Wiik et al. 1996; 


|Schrijver|2001||de Groof et al.|2005|[Muller et al.|2005a| b In or¬ 

der to obtain information about flows at other temperatures, we 
inspected co-spatial and co-temporal SDO/AIA data and found 
material flowing towards the sunspot in all channels. In the right 


panel of Fig. [d] we show the intensity in the 304 A passband, 
integrated across the path marked by the green line in the AIA 


304 A context image, as a function of the distance from point 
“C”. We find plane-of-sky inflows towards the sunspot reach¬ 
ing velocities of the order of the (vertical) downflow observed 
by IRIS (indicated by the continuous line of 90km s~'). One 
could speculate that these flows feed the downflow plume. How¬ 
ever, the signature of the inflowing material in the AIA channels 
disappears before reaching the light-bridge which is marked by 
the blue dashed lines in the time-distance plot and by the corre¬ 
sponding tick marks in the AIA 304 A context image in Fig. [4 
leaving the central part of the umbra dominated by the umbra 
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Fig. 7. IRIS 1400 A slit jaw image of the same area of AR 11836 as 
the HMI image in the lower left panel of Fig. [7] A time series of these 
images is available as a movie in the online edition. 


oscillations. These umbral oscillation are particularly strong in 
the 304 A, 171 A, 211 A passbands and seem to originate from 
the center of the umbra where the satellite feature is located, as 
already mentioned above. 

In addition, we inspected the IRIS 1400A slit-jaw data (see 
movie attached to Fig. |7J, which also shows flows towards the 
sunspot. Again, as seen in the AIA 304 A data, the flows become 
invisible at the outer boundary of the penumbra. No flows or 
loops can be detected that connect to the core of the umbra. 

Since none of these flows apparently connects to the center 
of the sunspot umbra where we observe the downflow in the tran¬ 
sition region, we cannot ascertain a link between the supersonic 
transition region downflow measured with the IRIS spectrograph 
and the inflows towards the spot observed with AIA and the IRIS 
slit-jaw cameras. Moreover, as in most observations of coronal 
rain, we could only track individual "blobs" falling along the 
magnetic field lines. These blobs, if they really fell through to 
the satellite region around point “C”, would appear as strongly 
intermittent events. The intermittency of coronal rain is also a 
key characteristic of models of coronal rain (M uller et al.| 2003, 
2004 2005b ). Further, coronal rain is often visible in chromo¬ 
spheric lines such as Her and Ca II, whereas the downflow we 
observe is visible only in transition region lines formed around 
10 5 K, with no trace of velocities in the cooler Cn and Mgii 
lines. We therefore believe that coronal rain is not consistent with 
most of the observational facts listed in Sec. [3] but we cannot ex¬ 
clude it with certainty. 

Another possible explanation, which we find more appeal¬ 
ing, in particular considering the remarkable stability of the su¬ 
personic downflow, is a siphon flow, along the lines discussed by 
e.g. |Cargill & Priest | ( |1980| ) and |Noci| | |1981) , who studied sta¬ 
tionary flow solutions in a coronal magnetic arch. In the follow¬ 
ing, we will refer to Noci’s paper since it makes predictions that 
are more relevant to our observations. His “class (vii)” solutions 


are of particular interest here. These solutions describe subsonic- 
supersonic steady flows with stationary shocks which adjust the 
flow to the boundary conditions in the second foot-point. The 
location and strength of the stationary shock is determined by 
the pressure ratio between the two foot-points. With decreasing 
pressure in the downflow foot-point the shock moves closer to¬ 
wards the foot-point with increasing jumps in temperature and 
density. We should point out, though, that Noci’s solutions are 
for coronal loops at coronal temperatures (> 1MK), whereas we 
observe the flow at transition region temperatures. However, we 
see no reason why these solutions should be limited to coronal 
temperatures and speculate that the supersonic downflow we see 
at transition region temperatures is the signature of a stationary 
shock in a siphon flow close to the foot-point in the sunspot um¬ 
bra. 

In this picture, the main component dominated by the um¬ 
bral 3-min shock waves is formed in the part of the transition 
region underneath the stationary shock of the siphon flow. As 
mentioned above, the strength of the stationary shock in Noci’s 
solutions depends on the pressure ratio between the foot-points. 
Due to the strong non-linear umbral 3-min oscillations, one 
cannot expect a constant pressure ratio between the two foot- 
points. One might therefore expect a modulation of the station¬ 
ary shock intensity by the underlying umbral 3-min shock waves. 
At minimum intensity in the main component, which coincides 
with minimum pressure and minimum temperature at the base 
of the downflow foot-point, Noci’s solutions would predict the 
strongest stationary shock in the siphon flow. This is exactly 
what we observe. As can be seen in the lower left panel of Fig.[6j 
the maximum intensity of the supersonic downflow coincides 
with minimum intensity of the main component. 

Without having identified the upflow foot-point, we are not 
in a position to conclusively claim the detection of a siphon 
flow at transition region temperatures; we can merely speculate. 
Along these lines it is interesting to note that about 21 hours 
after the observations presented here, the AIA 171 A channel 
shows prominent coronal loops connecting the dark quarter of 
the sunspot umbra directly to the trailing plage of opposite po¬ 
larity (Fig. [8]). Could there have been a similar loop at cooler, 
transition region temperatures during our observations? 

In order to get some further insight into the evolution and 
extent of the downflow, we inspected all IRIS observations of 
AR 11836 on 2015 September 2 and 3. In addition to the data 
set studied in detail in this paper, there are 7 more data sets of 
AR 11836: 5 large coarse 64-step rasters with 2" step width, and 
2 two-step sparse rasters. One of these data sets, the large coarse 
raster starting on 2013-09-02Tll:56:35 indeed shows clear sig¬ 
natures of a supersonic downflow for 4 steps (#35-38) during all 
5 raster scans. However, this downflow occurs at a different loca¬ 
tion in the sunspot, extending from the southern part of the most 
prominent light-bridge visible in Fig.[l](x ~ 4 Mm, y ~ -2 Mm) 
into the penumbra to the West. The measured downflow speed 
is smaller (« 75 kms -1 ). Since there are only 5 scans (with a 
scan cycle time of 188 s), it is difficult to infer information about 
temporal variations and its relation to the shock-wave dominated 
main component. All we can say is that the feature lasted for 
more than 12 min, has a similar extent along the slit as the down¬ 
flow studied here and extended over at least 8" in the direction 
perpendicular to the slit. The large coarse raster scan following 
the observations studied in this paper (starting at 18:29:35 UT) 
also shows some hints of a downflow, but with a much weaker 
signal and not in consecutive frames, as if it had split up. None of 
the other four data sets from September 2-3 shows any signs of 
a supersonic downflow. In a future paper, in which we will study 
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Fig. 8. AIA 171 A image of the same area as in Fig. [I] taken a day 
later on Sep 3, 2013 14:20 UT. Note the coronal loops (indicated by red 
arrows) connecting the dark quadrant of the spot directly to the trailing 
plage of opposite polarity. 


the range of main properties, locations and rate of occurrence of 
this phenomenon, we will analyse these data sets in more detail. 

From the measured density (cf. Section 3.3.2| > of the down¬ 
flow and its speed, v ~ 90 km s 1 , we can estimate the mass flux 
as pv — N e m p (Nu/N e ) v, where p is the mass density and m p 
is the proton mass, which yields « 5 x 10~ 7 gctrT 2 s _I assum¬ 
ing Nu/N e = 0.83. This is a substantial flux. It would evacuate 
the overlying corona on time scales of the order of 10 s. Even 
a prominence (column number density of 1 x 10 19 cm -2 ) would 
be drained in about 40 s. For the interpretation of the flow as a 
siphon flow, such a high mass flux does not pose a problem. 


5. Conclusions 


We have observed a small-scale, supersonic downflow of about 
90kms 1 in the transition region above the sunspot umbra of 
AR 11836. The downflow is not visible in the chromospheric 
lines, which only show an intensity enhancement at the loca¬ 
tion of the downflow. The downflow shows up as red-shifted, 
well-separated “satellite” lines of the Si iv and O iv transition re¬ 
gion lines and is remarkably steady over the observing period 
of nearly 80 min. The satellite lines do not participate in the 3- 
min shock wave Doppler maneuvres of the main component, but 
show a weak anti-correlated signature of the 3-min oscillations 
in intensity. We estimate the mass flux in the stationary flow to 
be « 5 x 10 7 gcrrT 2 s 1 . This mass flux is relevant for the atmo¬ 
spheric structure, as it would evacuate an overlying quiet coronal 
column in about 10 s. 

We interpret these findings as evidence of a stationary termi¬ 
nation shock of a supersonic siphon flow in a cool loop rooted 
in the central umbra of the spot, and we speculate that the super¬ 
sonic flows seen by |Dere|(|1982|),|Nicolas et aT7| ( | 1 982[ >, |Bryni IcD] 
sen et al. ( 2001| >, and Brynildsen et al. ( 2004[ ) are caused by the 
same phenomenon. We surmise that the 3-min shock waves vis¬ 
ible in the main component of the lines form in the transition 
region underneath the stationary shock and represent a perturba¬ 
tion to the boundary condition of the stationary flow. This would 
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explain the anti-correlation between the intensity of the main and 
the red-shifted satellite components. 

Models of these supersonic siphon flows with termination 
shocks were developed by Cargill & Priest| (JT980 ) and |Noci| 
(1981) for hot, coronal loops, although they should be appli¬ 
cable to transition region temperatures as well. There are nu¬ 
merical studies on siphon flows in coronal loops (e.g.: |Mariskaj 


1988 


Spadaro et al. |1991| Orlando et al.||1995a ), some includ¬ 


ing shocks (e.g Orlando et al.|1995b) , but we are not aware of 
any numerical studies reproducing shocks of the magnitude we 
are observing at transition region temperatures. Numerical stud¬ 
ies of flows in the transition region (e.g. Antiochosfl984 1 Peter 


|et al.|[2004] |2006[ |Hansteen et al .|[20 1 0[ > were usually targeted 
at reproducing the average red-shifts observed at those tempera¬ 
tures, and thus do not predict strong flows and shocks. It would 
therefore be interesting to consider numerical studies of solu¬ 
tions of the class studied by Noci, but for cool loops that never 
reach coronal temperatures. 

There are numerous observations in the photosphere, chro¬ 
mosphere and corona that have been interpreted as siphon flows 


(e.g. 

Rueedi et al. 1992 

[Spadaro et al. 2000 Uitenbroek et al. 

12006 

Beck et al. (2010; 

lethge et al. 2012). What are their rela- 


tions, if any, to the supersonic downflows measured in the transi¬ 
tion region over sunspot umbrae? Observations should also clar¬ 
ify whether steady-state supersonic downflows in the transition 
region are a common phenomenon, visible only in certain cir¬ 
cumstances, or whether they are a rather rare phenomenon. If 
they are a regularly occurring phenomenon, what is the range of 
their main properties? We find it intriguing, that the events with 
a single, clearly separated red-shifted component such as those 
observed by Dere ( 1982|), Nicolas et al.| ( |1982] >, one of the events 


of Brynildsen et a 


(200 


J, and this work, have approximately 
the same speed (~ 90kms 1 ), which is at the high end of the av¬ 
erage downflow speeds measured for 5 sunspots by Brynildsen 
|etafid2005) . 


Acknowledgements. IRIS is a NASA small explorer mission developed and op¬ 
erated by LMSAL with mission operations executed at NASA Ames Research 
center and major contributions to downlink communications funded by ESA and 
the Norwegian Space Center (NSC). CHIANTI is a collaborative project involv¬ 
ing George Mason University, the University of Michigan (USA) and the Uni¬ 
versity of Cambridge (UK). We gratefully acknowledge helpful discussions with 
Bart De Pontieu, Joe Gurman, Jim Klimchuk, Daniele Spadaro, Han Uitenbroek. 
This work has also benefited from discussions at the International Space Science 
Institute (ISSI) meeting on “Heating of the magnetized chromosphere” from 5- 
8 January, 2015, where many aspects of this paper were discussed with other 
colleagues. We also thank an anonymous referee for useful comments. 


References 

Antiochos, S. K. 1984, ApJ, 280, 416 
Avrett, E. H. & Loeser, R. 2008, ApJS, 175, 229 

Bartoe, J.-D. E, Brueckner, G. E., Nicolas, K. R., et al. 1979, MNRAS, 187,463 
Beck, C., Tritschler, A., & Woger, F. 2010, Astronomische Nachrichten, 331, 
574 

Bethge, C., Beck, C., Peter, H„ & Lagg, A. 2012, A&A, 537, A130 
Borrero, J. M. & Ichimoto, K. 2011, Living Reviews in Solar Physics, 8, 4 
Brekke, P, Brynildsen, N., Kjeldseth-Moe, O., Maltby, P., & Bartoe, J.-D. F. 

1991a, Advances in Space Research, 11, 251 
Brekke, P., Kjeldseth-Moe, O., Bartoe, J.-D. F., & Brueckner, G. E. 1991b, ApJS, 
75, 1337 

Brekke, P., Kjeldseth-Moe, O., & Brueckner, G. E. 1990, Ap&SS, 170, 135 
Brosius, J. W. & Landi, E. 2005, ApJ, 632, 1196 
Brynildsen, N., Maltby, P., Brekke, P„ et al. 1998, ApJ. 502, L85 
Brynildsen, N., Maltby, P., Kjeldseth-Moe, O., & Wilhelm, K. 2001, ApJ, 552, 
L77 

Brynildsen, N.. Maltby, P., Kjeldseth-Moe, O., & Wilhelm, K. 2004, ApJ, 612, 
1193 

Cargill. P. J. & Priest, E. R. 1980, Sol. Phys., 65, 251 










































































Straus et al.: A Steady-state Supersonic Downflow in the TR above a Sunspot Umbra 


de Groof, A., Bastiaensen, C., Muller, D. A. N., Berghmans, D., & Poedts, S. 
2005, A&A, 443, 319 

De Pontieu, B., Title, A. M., Lemen, J. R.. et al. 2014, Sol. Phys., 289, 2733 
Dere, K. P. 1982, Sol. Phys., 77, 77 

Dere, K. P, Landi, E., Mason, H. E„ Monsignori Fossi, B. C., & Young, P. R. 
1997, A&AS, 125. 149 

Doschek, G. A. & Feldman, U. 2004, ApJ, 600, 1061 

Dudlk, J., Del Zanna, G., Dzifcakova, E., Mason, H. E., & Golub. L. 2014, ApJ, 
780, L12 

Dwivedi, B. N. & Gupta, A. K. 1992, Sol. Phys., 138, 283 
Gurman, J. B. 1993, ApJ, 412, 865 
Hale, G. E. 1908, ApJ, 28,315 

Hansteen, V., De Pontieu, B., Carlsson, M., et al. 2014, Science, 346, 315 
Hansteen, V. H., Hara, H., De Pontieu, B., & Carlsson, M. 2010, ApJ, 718, 1070 
Jordan, C., Brueckner, G. E., Bartoe, J.-D. F., Sandlin, G. D., & Vanhoosier, 
M. E. 1978, ApJ, 226, 687 

Judge, P. G„ Carlsson, M„ & Stein, R. F. 2003, ApJ, 597, 1158 
Kjeldseth-Moe, O., Brynildsen, N., Brekke, P., et al. 1988, ApJ, 334, 1066 
Kjeldseth-Moe, O.. Brynildsen, N., Brekke, P, Maltby, P., & Brueckner, G. E. 
1993, Sol. Phys., 145, 257 

Kleint, L„ Antolin, P„ Tian, H„ et al. 2014, ApJ, 789, L42 
Landi, E., Young, P. R., Dere, K. P„ Del Zanna, G., & Mason, H. E. 2013, ApJ, 
763, 86 

Lemen, J. R., Title, A. M., Akin, D. J., et al. 2012, Sol. Phys., 275, 17 

Mariska, J. T. 1988, ApJ, 334, 489 

Mariska, J. T. 1992, Cambridge Astrophysics Series, 23 

Muller, D. A. N., de Groof, A., de Pontieu, B., & Hansteen, V. H. 2005a, 600, 30 
Muller, D. A. N„ De Groof, A., Hansteen, V. H„ & Peter, H. 2005b, A&A, 436, 
1067 

Muller, D. A. N„ Hansteen, V. H„ & Peter, H. 2003, A&A, 411, 605 
Muller, D. A. N„ Peter, H., & Hansteen, V. H. 2004, A&A, 424, 289 
Nicolas, K. R., Bartoe, J.-D. F., Brueckner, G. E., & Kjeldseth-Moe, O. 1982, 
Sol. Phys., 81,253 
Noci, G. 1981, Sol. Phys., 69, 63 
Orlando, S., Peres, G., & Serio, S. 1995a, A&A, 294, 861 
Orlando, S., Peres, G., & Serio, S. 1995b, A&A, 300, 549 
Pesnell, W. D., Thompson, B. J., & Chamberlin, P. C. 2012, Sol. Phys., 275, 3 
Peter, H„ Gudiksen, B. V., & Nordlund, A. 2004, ApJ, 617, L85 
Peter, H„ Gudiksen, B. V., & Nordlund, A. 2006, ApJ, 638, 1086 
Rempel, M. & Schlichenmaier, R. 2011, Living Reviews in Solar Physics, 8, 3 
Roussel-Dupre, R., Francis, M. H„ & Billings, D. E. 1979, MNRAS, 187, 9 
Rueedi, I., Solanki, S. K„ & Rabin, D. 1992^ A&A, 261, L21 
Scherrer, P. H., Schou, J.. Bush, R. I., et al. 2012, Sol. Phys., 275, 207 
Schrijver, C. J. 2001, Sol. Phys., 198, 325 

Spadaro, D., Antiochos, S. K., & Mariska, J. T. 1991, ApJ, 382, 338 
Spadaro, D., Lanzafame, A. C., Consoli, L., et al. 2000, A&A, 359, 716 
Tian, H., Curdt, W„ Teriaca, L.. Landi, E., & Marsch, E. 2009, A&A, 505, 307 
Tian, H., DeLuca, E., Reeves, K. K., et al. 2014a, ApJ, 786, 137 
Tian, H., DeLuca, E. E., Cranmer, S. R., et al. 2014b, Science, 346, A315 
Uitenbroek, H., Balasubramaniam, K. S., & Tritschler, A. 2006, ApJ. 645, 776 
Wiik, J. E., Schmieder, B., Heinzel, P., & Roudier, T. 1996, Sol. Phys., 166, 89 
Yurchyshyn, V., Abramenko, V., & Kilcik, A. 2015, ApJ, 798, 136 


Article number, page 11 of 1 111 


